arXiv:1507.04447vl [astro-ph.GA] 16Jul2015 


Mon. Not. R. Astron. Soc. 000 , ??-?? (2015) Printed 17 July 2015 (MN style file v2.2) 


Galaxy And Mass Assembly (GAMA): the effect of close 
interactions on star formation in galaxies 

L. J. M. Davies^*, A. S. G. Robotham^, S. P. Driverb^^ M. Alpaslan^, I. K. Baldry^, 

J. Bland-Hawthorn®, S. Brough^, M. J. I. Brown^, M. E. Cluver^®, M. J. Drinkwater^\ 
C. Foster®, M. W. Grootes^, I. S. Konstantopoulos®, M. A. Lara-Lopez^, 

A. R. L6pez-Sanchez®d2^ j, Loveday^^, M. J. Meyer^, A. J. Moffett^, P. Norberg^^, 

M. S. Owers®’^^, G. G. Popescu ^^^6^ P)g Propris^^, R. Sharp^^, R. J. Tuffs^, 

L. Wang^^’^^, S. M. Wilkins^^, L. Dunne^^’^®, N. Bourne^®, M. W. L. Smith^^ 

^ ICRAR, The University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia 
^ SUPA, School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews, Fife, KY16 OSS, UK 
^ NASA Ames Research Centre, N232, Moffett Field, Mountain View, CA 94035, United States of America 
Astrophysics Research Institute, Liverpool John Moores University, IC2, Liverpool Science Park, IfO Brownlow Hill, Liverpool, L3 5RF, UK 
^ Sydney Institute for Astronomy, School of Physics A28, University of Sydney, NSW 2006, Australia 
® Australian Astronomical Observatory, PC Box 915, North Ryde, NSW 1670, Australia 
^ Max Planck Institut fuer Kernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany 
^ School of Physics and Astronomy, Monash University, Clayton, Victoria 3800, Australia 

^ Instituto de Astronomia, Universidad Nacional Autonoma de Mexico, A.P. 70-264, 04510 Mexico, D.F., Mexico 
University of the Western Cape, Robert Sobukwe Road, Bellville, Cape Town 7535, South Africa 
School of Mathematics and Physics, University of Queensland, Brisbane, QLD 4072, Australia 
Department of Physics and Astronomy, Macquarie University, NSW 2109, Australia 

Astronomy Centre, Department of Physics and Astronomy, University of Sussex, Falmer, Brighton BNl 9QH 
Institute for Computational Cosmology, Department of Physics, Durham University, Durham, DHl 3LE, UK 
Jeremiah Horrocks Institute, University of Central Lancashire, PRl 2HE Preston, UK 
The Astronomical Institute of the Romanian Academy, Str. Cutitul de Argint 5, Bucharest, Romania 
Finnish Centre for Astronomy with ESC, University of Turku, Vaisdldntie 20, Piikkid, 21500, Finland 

Research School of Astronomy & Astrophysics, Mount Stromlo Observatory, Cotter Road, Weston Creek, ACT 2611, Australia 
Department of Physics and Astronomy, University of Canterbury, Private Bag 4800, Christchurch, 8 I 4 O, New Zealand 
Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ, UK 
School of Physics & Astronomy, Cardiff University, The Parade, Cardiff CF24 3AA, UK. 

SRON Netherlands Institute for Space Research, Landleven 12, 9747 AD, Groningen, The Netherlands. 


Accepted: May 2015 


ABSTRACT 

The modification of star formation (SF) in galaxy interactions is a complex process, 
with SF observed to be both enhanced in major mergers and suppressed in minor pair 
interactions. Such changes likely to arise on short timescales and be directly related 
to the galaxy-galaxy interaction time. Here we investigate the link between dynami¬ 
cal phase and direct measures of SF on different timescales for pair galaxies, targeting 
numerous star-format ion rate (SFR) indicators and comparing to pair separation, indi¬ 
vidual galaxy mass and pair mass ratio. We split our sample into the higher (primary) 
and lower (secondary) mass galaxies in each pair and find that SF is indeed enhanced 
in all primary galaxies but suppressed in secondaries of minor mergers. We find that 
changes in SF of primaries is consistent in both major and minor mergers, suggesting 
that SF in the more massive galaxy is agnostic to pair mass ratio. 

We also find that SF is enhanced/suppressed more strongly for short-time duration 
SFR indicators {e.g. Ha), highlighting recent changes to SF in these galaxies, which 
are likely to be induced by the interaction. We propose a scenario where the lower 
mass galaxy has its SF suppressed by gas heating or stripping, while the higher mass 
galaxy has its SF enhanced, potentially by tidal gas turbulence and shocks. This is 
consistent with the seemingly contradictory observations for both SF suppression and 
enhancement in close pairs. 
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1 INTRODUCTION 

Interactions are a key process in the evolution of galaxies in 
the Universe {e.g. White fc Frenk||1991 1. As galaxies interact 
and merge hierarchically over cosmic timescales, these inter¬ 
action events leave strong imprints on the galaxies involved, 


interaction. However, by contrast Li et al. (20081 also find 


modifying their morphology {e.g. Conselice et al. 2003 2005 


Lotz et al.||2008a|bl IMortlock et aL||2013 

, or see review in 

Gonselice (2014)), triggering AGN {e.g. 

Kazantzidis et al. 

2005| IHopkins et al.|2008| |Medling et al.|2013| however c.f. 

Villforth et al.|2014l, varying gas fractions {e.g. Kazantzidis 

et al.|2005||Lotz et al.|2010||Ueda et al.|2012||2013P and sig- 

nificantly altering their star-formation (SF) history {e.g. 

Di 

Matteo et al.|20071|Shi et al.|20091 

Wong et al. 120111 |Bour- 

naud et al.|2011| Patton et al.|2013 

[Robotham et al.|2013|. 


Such merger events are ubiquitous throughout the Universe 
and as such, probing the details of the merger process is 
key to our understanding of galaxy formation and evolution 
{e.g. Casteels et al.|26l4i. For example, star-formation (SF) 


initiated by interactions could well be a key factor in trans¬ 
forming blue star-forming discs into red, passively evolving 
quiescent spheroids. 

However, the details of how SF responds to galaxy in¬ 
teractions are somewhat vague and are likely to be strongly 
linked to pair mass ratio and pair separation. Tradition¬ 
ally, observational studies of close pairs have highlighted 
strong evidence that SF is enhanced through interactions 
and that the strongest enhancement occurs in the closest 


pairs (those at projected separations of <30kpc, e.g. Elli- 
son et aL||2008 Freedman Woods et al.|[2010|. In addition. 


detailed multi-wavelength analyses of blue compact dwarf 
galaxies have shown that, in the majority of the cases, in¬ 
teractions with or between low-luminosity dwarf galaxies or 
HI clouds are the triggering mechanism of their strong star¬ 


formation activity (Lopez-Sanchez 2010| 

Lopez-Sanchez et 

al. 20121. More recently, Patton et al.| ( 

20131) have stud- 


ted a large sample of star-forming galaxies in pairs taken 
from the Sloan Digital Sky Survey (SDSS) and find clear 
evidence for enhanced, emission-line derived SF out to sig¬ 
nificant pair separations (~ 150 kpc). They also find that 
the enhancement in SF is inversely proportional to pair sep¬ 
aration, with the closest pairs displaying the largest en¬ 
hancement. This result is echoed by |Wong et al.| ( |2011[ ) 
who perform a similar analysis for isolated pair galaxies 
in the PRIsm MUlti-object Survey (PRIMUS) but deter¬ 
mine star-formation rates (SFRs) using attenuation cor¬ 
rected Ultraviolet(UV)-optical colours. They find pairs at 
separations of < 50 kpc show bluer far-ultraviolet minus r- 
band (FUV-r) colours than a control sample of non-pair 
galaxies (indicative of larger SFRs). Moreover, this SFR en¬ 
hancement is more pronounced at closest pair separations 
(< 30kpc). [Scudder et al.| ( |2012[ ) perform a similar analy¬ 
sis on SDSS pairs and, in addition to SF enhancement, find 
that SF is most strongly enhanced in major merger systems, 
hinting that pair mass ratio is significant in the modification 
of SF history through galaxy interactions. In a distinct but 
complementary approach, [Owers et al.| ( |2007[ ) find that 2df 
Galaxy Redshift Survey (2dfGRS) selected star-busts are 
more likely to have a neighbour of comparable brightness 
within 20 kpc, and 30% of star-bust galaxies show morpho¬ 
logical signatures of interactions and mergers - suggesting 
that the SF in these systems is induced by a galaxy-galaxy 


strong enhancement of SF in interactions but suggest that 
there is little correlation between this enhancement and the 
relative luminosity of the interacting galaxies (and therefore, 
pair mass ratio). They do also find that the enhancement of 
SF is stronger for lower mass systems. By splitting individ¬ 
ual interacting galaxies into high (< log[M*/M q\ >= 10.6) 
and low (< log[M*/M q] >= 9.72) stellar mass samples, 
they find that at a given pair separation, SF is enhanced 
more strongly in low mass than high mass pair galaxies. 

Numerical simulations offer further insight into the 
complexities of these interaction SF processes. |Di Matteo| 
|et al.| ( |2007| model several hundred galaxy close pair in¬ 
teractions (< 20 kpc) for various morphological classes, and 
find that while SF is primarily enhanced in galaxy interac¬ 
tions, mergers do not always trigger starbursts, and galaxy 
interactions are not always sufficient to convert high gas 
masses into new stars. They also highlight that the amount 
by which SF is enhanced is anticorrelated with pair sepa¬ 
ration on small scales {i.e. galaxies which have a very close 
passage produce the lowest bursts of SF), as well as the am¬ 
plitude of tidal forces {i.e. pairs that undergo less intense 
tidal forces can preserve higher gas masses for future SF 
during the merger). These simulations suggest that SF en¬ 
hancement is not ubiquitous in interactions and that SF is 
likely to vary as a function of pair separation and pair mass 
ratio - higher mass ratios will induce more significant tidal 
effects in the lower mass galaxy in the pair removing the 
bulk of its gas and as such, starving SF. 

Recently we have seen tentative observational evidence 
in support of these simulation predictions for close pairs. 
Robotham et al. (20131 studied L* galaxies in closely in¬ 


teracting pairs taken from the Galaxy And Mass Assembly 
(GAMA) Galaxy Group Catalogue (G^C), and investigate 
variations in SF as a function of pair mass ratio. They find 
that primary (higher mass) pair galaxies have their SF en¬ 
hanced, but only when the pair mass ratio is close to 1. 
While the secondary (lower mass) galaxy has its SF sup¬ 
pressed relative to equivalent stellar mass non-pair galax¬ 
ies. In addition, De Propris et al. (20141 find that close 


pairs galaxies are consistently redder than a similarly se¬ 
lected comparison sample, and suggest that these systems 
have been “harassed” in multiple previous passes prior to the 
current close interaction. These results give the first tenta¬ 
tive observational evidence that SF in interacting systems is 
more complicated than the traditional evidence would sug¬ 
gest, and highlights that pair mass ratio may play a vital 
role in the changes to SF induced by an interaction. 

However, questions still remain regarding the details of 
SF in interacting systems, such as how does SF proceed in 
interactions as a function of pair separation when consider¬ 
ing pair samples split by merger ratio, primary or secondary 
status and individual galaxy mass? A potentially more sig¬ 
nificant question is over what timescales do variations in the 
SFR in interacting systems become apparent? If we measure 
SFRs using observational tracers which probe different SF 
timescales, will we see differences in observed SFRs? SF en¬ 
hancement/suppression in interactions is likely to occur on 
short timescales, and hence, may only be apparent in SFR 
measures which probe short periods of the galaxy’s SF his¬ 
tory. 

In this paper, we further investigate SF in closely in- 
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teracting pairs taken from the G^C catalogue. We split our 
sample on merger ratio, primary or secondary status and 
galaxy mass, and investigate SF characteristics as a function 
of pair separation. We determine SFRs for each galaxy using 
multiple nebular emission line and continuum fluxes, which 
probe the galaxy’s SF history over different timescales, and 
compare these methods in order to identify short timescale 
variations in SF, once again as a function of pair separation. 
In this manner we identify factors driving the modification 
of star-formation rates in interacting galaxies. In Section 
we outline the GAMA survey and pair catalogue used in 
this work, Section discusses the SFR indicators used, in 
Section we highlight the effects of interactions on SF in 
galaxies and in Section we summarise our results. 

Throughout this paper we use a standard ACDM cos¬ 
mology with Ho = 70kms“^ Mpc“^, Ha = 0.7 and Hm = 0.3. 


2 DATA 

2.1 GAMA and the Pair Catalogue 


The GAMA survey is a highly complete multi-wavelength 
database (Driver et al. poTT|) and galaxy redshift {z) survey 


I Baldry et al.||2010[ Hopkins et al.|[2013| [Liske et al.||2014 1 
covering 280degro a main survey limit of tab < 19.8 mag 
in three equatorial (GOO, G12 and G15) and two southern 
(G02 and G23) regions. The spectroscopic survey was un¬ 


dertaken using the AAOmega fibre-fed spectrograph (Sharp 
|et al.|[^006| [Saunders et al.||2004[) in conjunction with the 


Two-degree Field (2dF, Lewis et al. 20021 positioner on 


the Anglo-Australian telescope and obtained redshifts for 
~250,000 targets covering 0 < a < 0.5 with a median red- 
shift of z ~ 0.2, and highly uniform spatial completeness 


(Baldry et al. 2010 Robotham et al. 2010 Driver et al. 
2011D. Full details of the GAMA survey can be found in 


Driver et al. (20111 and Liske et al. (20141. In this work 


we utilise the first 5 years of data obtained and frozen for 
internal team use, referred to as GAMA 11. 

Firstly, to minimise Active Galactic Nuclei (AGN) con¬ 
tamination in our sample, which may potentially bias our 
SFR estimates, we exclude optically bright AGN using the 


BPT diagnostic diagram (Baldwin et al. 19811. We select 


sources from the GAMA II spectral line catalogue which 
have all BPT diagnostic lines detected at S/N > 3, and have 
emission-lines which do not lie at the edge of the spectral 
range (and may have poor line measurements). We then ap¬ 
ply the BPT AGN-bcomposite source selection line of jKauff^ 
mann et al. (20031: 


logio([OIII]/H^) > 


0.61 


-b 1.3, 


( 1 ) 


logio([NII]/Ha)-0.05 

to identify galaxies which may have their SFR contaminated 
by AGN (Fig. [^. Galaxies which meet the BPT AGN se¬ 
lection are coloured green and those which are potentially 
composite AGN-bSF sources are coloured orange. We find a 
total of 2,486 galaxies from the full GAMA H sample which 
potentially contain an optically bright AGN. 

In order to select systems which are undergoing a close 
interaction, we utilise the GAMA G^C catalogue which in¬ 
cludes the identification of all galaxy pairs ([Robotham et al. 


2011 also see Robotham et al. 2012 2013 2014 b Briefly, 



Figure 1. The BPT diagram for all GAMA II sources with all 
diagnostic emission lines detected at S/N > 3. The AGN dividing 
line of [Kewley et ah] | [2013[ l is displayed as the black dashed line 
while the composite dividing line of [Kauffmann et ah] | [2003[ l is 
shown as the black solid line. Sources which potentially contain 
AGN are displayed as green or orange. 


pairs are selected on physical projected separation, rsep-proj 
(for the cosmology given in Section and radial velocity 
separation, Ugep-rad. In this paper we use systems which 
meet the pair criteria of Robotham et al.[([2012| 2013 2014 b 


-Prl 00 v 500 — rgep —proj ^ 100 kpC A Ugep —rad ^ 500 km S 

This sample consists of a total of 37,679 galaxies in 
pairs. We further restrict our sample to pairs at z < 0.3 in 
order to minimise the impact of redshift evolution in our 
galaxies and biases towards the identification of high mass 
pairs at the high redshift end, giving a total of 33,832 pair 
galaxies. We then match all pair galaxies to the GAMA H 
panchromatic photometry catalogue (Driver et al in prep) 
and use rest-frame photometry derived from a refactored 


implementation of the InterRest algorithm ( Rudnick et al. 
2003 Taylor et al.[[200'9 l, coupled with the empirical set of 
galaxy template spectra of Brown et al. (2014 b We exclude 
any galaxy classed as an AGN using the classifiers described 
above, leaving 32,468 pair galaxies. 

We use the stellar masses derived from the ugriZJH 


photometry for all GAMA II galaxies (Taylor et al. 20111 


to assign classes to each individual galaxy. Firstly, we de¬ 
termine whether each galaxy is the primary or secondary 
system in the interacting pair (by mass). Secondly, we cal¬ 
culate the pair mass ratio and identify each system as being 
either part of a major merger (mass ratio <3:1) or minor 
merger (mass ratio >3:1). Note that we also split our minor 
merger sample further into 3:l<mass ratio<6:l and mass 
ratio>6:l subsamples, but see little difference between each 
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of these populations. As such, we do not discuss them in this 
work. 

A potential caveat to the selection above, is that we 
may be biased towards pair systems where both galaxies are 
star-forming, specifically at the low mass and/or high red- 
shift end - where we are not mass complete. For example, by 
purely applying a redshift selection we are more sensitive to 
higher stellar mass and, potentially star-forming, galaxies. 
As such, low mass systems will only be identified as being 
in an interaction if both pair galaxies are high mass and/or 
highly star-forming. Likewise, if we were to only apply a 
mass selection over the full sample, we do not take into 
account the broad redshift range of our selection and any 
evolution of the star-forming properties of galaxies across 
this extensive look-back time. As such, our pair sample will 
be biased to wards specihc galaxy populations and we will 
dilute any observed trends by not taking into account the 
global change in star-formation properties of galaxies with 
redshift. In order to remove this affect, we compare all of our 
pair galaxy samples to a mass and redshift matched control 
sample. The control sample contains the same mass bias and 
redshift evolution as the pairs sample, and as such removes 
any dependancy on this effect (GAMA is essentially 100% 
complete in the r-band, thus sources will also not drop out 
of our pairs selection preferentially to our control sample 
due to the effect of the interaction on SF - where the r-band 
will not be strongly affected). As such, we can directly com¬ 
pare the effects of an interaction to non-interacting galaxies 
at the same stellar mass and redshift. Our control sample is 
defined further in Sec. |4.2[ In this paper we only display our 
results in comparison to this control sample, and as such are 
not biased by our mass incompleteness. 

We note that in defining our control sample based on 
stellar mass and redshift, we are not taking into account the 
different galaxy clustering statistics of passive and active 
systems. Passive galaxies are more strongly clustered (e.g. 


Zehavi et al. 20111 and as such, more likely to be in pairs 
(see |Lin et al.|2008| |de Ravel et al.|2009[ for the fraction of 
quiescent galaxies in major mergers at z ~ 0). This could 
potentially bias our control sample towards a higher fraction 
of actively star-forming systems. However, we do not wish 
to perform any sample selection based on star-formation di¬ 
agnostics, as the increased fraction of passive systems in 
pair galaxies may in fact be caused by the galaxy interac¬ 
tions - the very the effect we wish to measure here. As we 
will see, the majority of sub-samples in our analysis show 
enhanced star-formation in close pairs, a result that would 
only be reduced by the potential biases discussed above. 
In addition, we see distinct differences between the primary 
and secondary galaxies of minor mergers in our sample. Any 
bias produced by different clustering statistics would affect 
both populations in the same manner. Hence, such a bias is 
unlikely to be driving this result. 

In all of the subsequent analysis in this paper, we con¬ 
sider each pair galaxy member individually, but split our 
samples into primary and secondary galaxies, and major and 
minor mergers based on the classihcation within their pair 
system. We also further split our samples by the stellar mass 
of each individual galaxy as noted above. For reference, Ta¬ 
ble shows the number of galaxies in each mass bin as a 
function of primary or secondary status and major or minor 
merger classification. 


Table 1. Number of pair sources (each pair galaxy is treated in¬ 
dividually) in each mass bin as a function of primary or secondary 
system status and major or minor merger classification. 


Mass Range 
(logic [Mq]) 

Total 

Primary 

Secondary 

Major 

Minor 

9.5-10.0 

4,606 

1,407 

3,199 

2,088 

2,517 

10.0-10.5 

8,795 

3,958 

4,837 

5,121 

3,670 

10.5-11.0 

10,910 

7,282 

3,628 

6,870 

4,040 

All 

24,311 

12,647 

11,664 

14,079 

10,227 


Short Duration 

<- 


Long Duration 


MagPhys O.tGyr 


H-alpha 



10 100 
Probes SF, Myr 


Figure 2. SF timescale probed by each SFR indicator used in 
this work, for normal disk galaxies. Sources of the emission are 
noted on the right hand side of the figure. We split our measures 
into short (< 100 Myr) and long (> 100 Myr) duration. See text 
for details of each SFR indicator. 


3 SF MEASURES OVER DIFFERENT 
TIMESCALES 

In this paper we consider SFRs determined via different ob¬ 
servables in order to probe SF on different timescales and 
equate them to recent variations in the galaxy’s instanta¬ 
neous SFR. Fig. shows a cartoon pictorial representa¬ 
tion of the timescales over which each SFR indicator (de¬ 
scribed in more detail below) probes in normal disk galaxies, 
and the source of the emission in each case. This figure is 
only intended to highlight the key differences in timescales 
over which each indicator probes, as the complex physics 
of each process is likely to place large variation on the true 
timescale of the emission arising from each source. As such, 
we do not constrain our SFR indicators further than split¬ 
ting them into long duration - probing timescales >100 Myr, 
far-infrared (FIR), UV-|-total infrared (UV-I-TIR) and mid- 
infrared (MIR), and short duration - SFR indicators probing 
only < 100 Myr, FUV, near-ultraviolet (NUV ), MAGPHYS 
(0.1 Gyr) and Ha. See Gilbank et al.| (2010| for a compar¬ 
ison of different SFR indicators in the SDSS stripe 82 and 
details of each SFR indicator used in this work below. 
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3.1 Long duration MIR/FIR Continuum SFRs 


Infrared (IR) emission at l-lOOO/rm arising from ‘normal’ 
star-forming galaxies is produced from three sources: photo¬ 
spheres and circumstellar envelopes of old stars undergoing 


mass loss {e.g. Melbourne et al. 20121, interstellar gas and 
dust heated by either bright OB stars in star-forming regions 
(warm dust) or the general stellar radiation field through¬ 
out the Inter Stellar Medium (ISM, cool dust - ‘cirrus’) - 
for example see the review by [Sauvage, Tuffs, fc Popescu 
(20051, Popescu et al. (20001 or more recently Xilouris et al. 


(20121 and references therein. Stellar sources of IR emission 


dominate at short wavelengths < 3pm and interstellar gas 
emission makes up just a few percent of the total IR output 
of galaxies. At 3— 100pm the bulk of the emission arises from 
warm dust locally heated by UV emission from young stars 
in star-forming regions, with cirrus emission dominating at 
> 60pm. As such, probing IR flux from star-forming galaxies 
in the 3 — 100pm range gives a reliable estimate for the ongo¬ 
ing SF {e.g. Calzetti et al. 2007| - the amount of flux emitted 
in the IR is directly related to the UV emission from newly 
formed stars). However, IR emission from dust requires sig¬ 
nificantly long timescales to become apparent and subsides 
slowly when SF is suppressed (see Kennicutt|1998 K As such, 
SFRs derived from the IR continuum probe SF over large 
timescales, giving an estimate of a galaxy’s SF history on 
timescales of > 100 Myr. 

In this paper we consider two IR continuum measures 
of SFR. Firstly, we determine the FIR SFRs derived from 
the 100pm flux provided for all GAMA sources as part of 
the Herschel Astrophysical Terahertz Large Area Survey 


(H-ATLAS Bales et al. 20101 and outlined in the GAMA 
II panchromatic data release (Driver et al in prep). =The 
100pm data comes from the H-ATLAS Phase 1 Data Re¬ 
lease (Valiante et ah, in prep), which provides PACS maps at 
100 and 160pm reduced using Scanamorphos. The GAMA- 
II photometry in the Herschel maps exploits an algorithm 


developed by Bourne et al. (20121 to optimally capture 


extended flux while accounting for blending in the low- 
resolution images. This algorithm consists of convolving the 
100pm map with a kernel given by the GAMA r-band- 
defined aperture smoothed with the 9” (FWHM) PSF of 
the PACS 100pm band. These smoothed apertures are cre¬ 
ated for all GAMA galaxies, and any that overlap in the 
map are down-weighted so that the 100pm flux in such pix¬ 
els is shared between the overlapping apertures. Fluxes in 
all other H-ATLAS bands are used for the UV-fTIR and 
Magphys SFRs (see below) and are produced in a similar 
manner, but using different PSFs - which are specific to the 


band in question. For more details see Bourne et al. (20121 


and Driver et al (in prep). We convert 100pm fluxes to SFRs 
using the tight correlation derived in |Davies et al.| ( [2014p for 
Virgo cluster galaxies: 


logio SFRmir(Mo yr ^) = 0.841ogioi^L22Mm(Lo) - 7.3. 

_ ^ 

We note that the correlations derived in ICluver et al.l 


(20141 compare WISE fluxes with Ha derived SFRs and as 


such, may be biased towards emission line-derived SFRs. 
However, we rescale all SFRs derived in the paper to have 
the same slope and normalisation as the SFRf/_r in order 
to compare each measure directly (see below). 


3.2 Short duration attenuation corrected UV 
Continuum and Magphys SFRs 

In contrast with the FIR and MIR, UV continuum derived 
SFRs probe a shorter timescale than those in the IR (< 
100Myr). UV continuum arises from hot, massive (M* > 
3Mq) O and B stars, and as such is a good tracer of more 
recent SF in galaxies (e.g. Kennicutt fc Evans||2012| ), with 
luminosity-weighted mean ages for a constant SFR predicted 
to be ~28 Myr for the GALEX FUV band, and ~80 Myr for 
GALEX NUV (Grootes et al. in prep). 

Hence, by comparing IR continuum and UV continuum 
SFRs we may be able to disentangle recent SF changes in a 
system - potentially due to interactions. 

However, direct estimates of the SFR from UV contin¬ 
uum luminosities are problematic as the UV emission arising 
from a galaxy is extremely sensitive to dust attenuation {e.g. 
Wang fc Heckman||1996 1. Previous studies have attempted 
to overcome this issue through estimating attenuation cor¬ 
rected SFRs by applying attenuation estimates based on the 
UV spectral slope (/3) and luminosity corrections, such as 
those derived by |Meurer et al. (19991 - e.g. Wijesinghe et al. 
(20111. Such corrections apply general scaling to all galaxies 


and may not be appropriate for specific galaxy classes {e.g. 
Wilkins et al.|20T2 l, severely biasing any results derived by 
such an analysis. 

Here we take a different approach and use the full Spec¬ 


tral Energy Distribution (SED) Magphys (da Cunha et al. 
20081 fits to the GAMA H galaxies (Driver et al in prep). In 


this upcoming work. Driver et al obtain the best SED fit to 
the full 21-band photometric data available to all GAMA H 
sources, simultaneously fitting the UV through FIR flux and 
obtaining the best-fit un-attenuated model spectrum. For 
our attenuation corrected UV luminosities we use the unat¬ 
tenuated UV emission from the best-ht Magphys model tem¬ 
plate convolved with the GALEX FUV/NUV filter curves. 
In this manner we use each galaxy’s full SED to estimate the 
attenuation correction to be applied to the UV luminosity 
for each individual galaxy instead of applying a general /3 
correction. We convert UV luminosity to at SFR using the 
calibrations given in |Wijesinghe et ah] (|2011[) : 


logic SFRFiR(MQyr ^) = 0.73logioLiooi.m(WHz ^) - 17.1. 

( 2 ) 

Secondly, we derive MIR SFRs using the Wide-held In¬ 
frared Explorer {WISE) data from the matched GAMA- 
WISE catalogue outlined in Cluver et al. (20141. We use 


WISE 22fj,m (W4) band huxes, which are not strongly af¬ 
fected by PAH emission, and the best-ht SFR correlation 


obtained in Cluver et al. (20141: 


SFRNuv(M0yr-i) = Lnuv(WHz i) 


SFRFuv(M0yr-i) = 


1.56 X 1021 


Lfuv(WHz-^ 
1.64 X 1021 ' 


( 4 ) 


( 5 ) 


The Magphys code also provides an estimate of the 
galaxy SFR averaged over the last 100 Myr using a best 
ht model for both the FIR and UV emission for a com- 
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bined estimate of both obscured and unobscured SF (here¬ 
after SFRo.i Gyr)- Hence, this provides an additional short 
timescale SFR estimate with which to compare to our long 
duration SFKfir- 


3.3 Long duration UV Continuum + Total 
Infrared Luminosity SFR 

We also use the combination of UV and total IR (TIR) lu¬ 
minosities as a star formation rate proxy broadly probing 
the last ~ 300 Myr of the galaxies star-formation history. 
As discussed above, UV emission arises directly from star¬ 
forming regions and probes SF on short timescales, while 
some fraction of this emission is absorbed and reprocessed 
by dust, being re-emitted in the FIR on longer timescales. 
As such, using a star-formation rate indicator which probes 
both the UV and FIR emission with give a relatively sta¬ 
ble, but broad timescale, measure of star-formation in our 
sample galaxies. We therefore sum both UV and TIR lumi¬ 
nosities to obtain a total star formation rate estimate, based 
on the bolometric luminosity of OB stars. 

We use the method outlined many high redshift studies 
(e.g.|Bell et al.|2005||Papovich et al.|2007||Barro et al.|2011|) 

of: 


SFRuv-i-tir(M0 yr ^) = 1.09 x 10 ^°[Lir -|-2.2Luv](L0). 

( 6 ) 

This prescription is the Bell et al. ( 2005|) recalibration 


of the relation from Kennicutt (19981, scaled for a Chabrier 


(20031 stellar IMF. Here, Lir is the total IR luminosity, 
integrated between 8-1000 fim. These values have been esti¬ 
mated by fitting the Chary & Elbaz (20011 galaxy templates 
to the WISE and Herschel photometric points. The bolomet¬ 
ric UV luminosity between 1216-3000 A, Luv, is estimated 
as 1.5:^/iy,2800, where fv ,2800 is the rest-frame luminosity at 
2800 A. For further details of this process see Taylor et al. 
(in prep). 


3.4 Short duration nebular Ha SFRs 


Hq photons arise from gas ionised by the stellar radiation 
held, and only stars with ages <20 Myr can contribute signif¬ 
icantly to this ionizing hux. Thus, Ha provides a direct mea¬ 
sure of the current SFR in galaxies (<10-20 Myr) which is 
largely independent of SF history (e.g. see Kennicutt|1998 1. 

For SFRhq we use emission line data from the GAMA H 
spectroscopic campaign (Owers et al. in prep), where aper¬ 
ture, obscuration and stellar absorption corrected Ha lumi¬ 
nosities are given by: 


SFRHc(M0/yr) = 


Lhc(WHz-^) 
1.27 X 1034 


( 8 ) 


One caveat to using Ha SFRs is that the aperture 
based spectroscopy only probes the central regions of nearby 
galaxies. However, in this paper we only investigate star- 
formation rates in comparison to a control sample of galaxies 
matched on mass and redshift (see below) and as such, both 
our pair and control sample should suffer the same aperture 
bias. 


We note that the GAMA stellar masses, UV-I-TIR and 
Magphys-based SFRs are calculated assuming a Ghabrier 
IMF, while our Ha, MIR, and FIR SFRs are calculated us¬ 
ing a Salpeter IMF, as such we scale all Salpeter IMF SFRs 


by a factor of 1.5 to account for this discrepancy (Dave 2008 
Driver et al.||20f3|. In addition, all results in this work are 


displayed relative to a control sample with SFRs derived in 
an identical manner. As such, differences been IMF assump¬ 
tions in in different indicators will not affect our results. 


3.5 Calibrating continuum SFRs 

As noted above, we recalibrate all continuum SFR measures 
to have the same slope and normalisation as SFRfjr. We 
apply this correction as all SFR indicators are derived us¬ 
ing vastly different calibrations, and different samples. If 
SFR indicators consistently derived SFRs in the same man¬ 
ner, we would find a similar slope and normalisation when 
comparing indicators. In Fig. we display our SFR indi¬ 
cators in comparison to SFRf/r for all sources detected 
at S/N > 2 in both the Herschel lOO/rm and WISE W4 
observations (the limiting bands in our selection). Glearly, 
not all indicators show the same slope and normalisation. 
These offsets are likely to be produced by the SF cali¬ 
bration, not the true galaxy population. As such, we ht 
the slope and normalisation of each distribution using the 
[R] multi-dimensional MGMG fitting procedure, Hyperfit 
(R obotham & Obreschkow, 2015) and recalibrate each SFR 
measure using the best fit parameter as follows: 


logic SFRcai = a X logic SFRoUg + b (9) 


where. 


SFRmir, a = 0.824 A b = 0.223 
SFRuv-ttir, a = 0.755 A b = 0.110 


Lhc = (EWh. -f EWc) X io-°-4(m.-34.i) 

3xl0i« /Fn^/Fn/iy-®® 

^ (6564.1(1-bz)2) 2.86 ) 


(7) 


and EWhq denotes the Ha equivalent width, EWc is 
the equivalent width correction for stellar absorption (2.5A 
for GAMA, [Hopkins et al.|[2013p , Mr is the galaxy r-band 
magnitude and Fua/FH/S is the Balmer decrement (see Gu- 
nawardhana et al.||201l1 for further details). Using this, 
SFRhci can be determined from jifennicuttj ( |1998| ), assum¬ 
ing a Salpeter IMF: 


SFRfuv, a = 0.610 A b = 0.491 
SFRnuv, a = 0.649 A b = 0.462 


SFRc.i Gyr, a = 0.595 A b = 0.562 
SFRhci , a = 1 a b = 0 

Note that we do not apply a correction to the Ha emis¬ 
sion line-derived SFR as there is large scatter in the distri¬ 
bution but the locus follows a 1:1 correlation. In fact, this 
scatter could be due to the very timescale variations we hope 
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SFR-MIR SFR-UV_TIR SFR-FUV SFR-NUV SFR-0.1Gyr SFR-Ha 


Figure 3. Comparisons between various SFR measures used in this work. The 1:1 relation is shown by the solid black line. Points are 
colour coded by source redshift, z. Hyperfit best fit lines to each correlation are displayed as the green lines with ztlcr errors on the 
normalisation. These fits are used to scale each distribution to the same slope and normalisation using equation]^ Note we do not scale 
SFRhci as the correlation between SFRs is poor. 


to identify in this work - short time scale variations will be 
most apparent throngh a comparison of Ha and FIR derived 
SFRs. We explore this further in Section [5.1| While we ap¬ 
ply these calibrations to form the main comparison in this 
work, we note that our conclusions hold true if no correc¬ 
tion is applied. A comparison of SFR indicators for the full 
GAMA II sample will be the subject of an upcoming paper 
(Davies et al in prep). 

One potential caveat to these scalings is that they are 
defined only for FIR detected sources, and may not be ap¬ 
plicable to the general galaxy population. To highlight this 
this is unlikely to cause significant bias in our results and 
that these scalings are appropriate. Fig. shows the main 
sequence of star-forming galaxies for a completely indepen¬ 
dent SFR indicator not used in this work. Here we use the 
extinction corrected NUV SFR (using the direct NUV fluxes 
and UV-spectral slope extinction correction, unlike the Mag- 
phys NUV fits used elsewhere in this work), which does not 
use any information from the FIR (this SFR indicator will 
be defined further in Davies et al in prep). For all GAMA H 
galaxies in our redshift and mass limited sample, we display 
sources with FIR detections (> 2cr) as blue points, and those 
without FIR detections as red points. Glearly, the main se¬ 
quence is consistent between FIR detected and undetected 
sources, suggesting that any correlations derived from the 
FIR detected sample is applicable to all galaxies. 

In the following analysis, unless otherwise stated, we use 
specific SFRs (sSFR) calculated by dividing all measured 
SFRs by stellar mass. This removes any mass dependance 
on observed SFR correlations. 


4 SF IN PAIR GALAXIES AS A FUNCTION 
OF PAIR SEPARATION 

In attempting to pin down the effect of close interactions on 
galaxy SFRs we have two measurements which can both be 
considered as a proxy for the stage of interaction which we 
are witnessing: i) pair separation and ii) SFRs on different 
timescales. In the former, to first order, we can equate the 
pair separation to stage of the interaction - in that distant 
pairs are more likely to be at an early stage and close pairs 
at a later stage of the interaction (we will discuss the caveats 
to this assumption in Section]^. Therefore, we may expect 
galaxies at large pair separations to display SFRs which are 



Stellar Mass, Msun 

Figure 4. The main sequence of star-forming galaxies for ex¬ 
tinction corrected NUV SFRs (see text for details). FIR detected 
sources are shown as blue points, while FIR undetected sources 
are shown as red points. Both populations are consistent and as 
such, correlations derived for the FIR detected sources are appli¬ 
cable to the full sample. 


not strongly affected by the interaction and close pairs to 
show the strongest effects. In the latter, we are directly mea¬ 
suring SFRs at different stages of the interaction. As such, 
we may expect long duration SFR indicators to probe an 
epoch prior to the modification of SF, while short duration 
indicators should once again show the strongest affect. 

In the following three subsections we separate out the 
variations with pair separation alone, considering each SFR 
indicator independently, and do not discuss the use of multi¬ 
ple indicators to highlight short duration timescale changes. 
However, when reading the following analysis we encourage 
the reader to keep in mind that pair separation is a proxy 
for timescale within the interaction. The following sections 
discuss the fine details of SFR variation in pairs and we di¬ 
rect the casual reader to Section [4.2.1 [ for a summary of the 
key results found in the following Section, and Section for 
a full summary of all of our results. 
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4.1 Comparisons between pair samples in 
Robotham et al. (2014) 

Initially, we compare different SFR indicators for our pair 
sample as a whole (not split on mass, pair mass ratio, etc) in 
order to highlight discrepancies in the results derived from 
previous studies probing SF in interacting systems. We split 
our pair systems on the three pair selection criteria outlined 


in Robotham et al. (20141: 


Pr20v500 = l"s 


< 20h ^kpc A Vs 


< 500 km s 


( 10 ) 


Pr50v500 = Vs 


< 50h ^kpc A Vs 


< 500 km s 


Rr- 100 t ;1000 = Vssp < 100 ^kpc A Usep < 1000 km s ^ 

with each sample excluding sources in the inner contained 
samples {i.e. galaxies in the Pr 20 i ;500 sample are not con¬ 
tained in the Prsousoo or Prioouiooo sample). Hereafter we 
define these samples as close-, intermediate-, and far-pairs 
respectively. We calculate the median logic [sSFR] for each 
indicator and normalise to the value for far-pairs - to high¬ 
light variations in sSFR between closely interacting and dis¬ 
tant systems. Fig. shows the difference between pair sam¬ 
ples for each SFR indicator, points are offset for each pair 
sample in the x-direction for clarity. Error bars display the 
standard error on the median for each sample, where for 
many points this is smaller than the plotted symbol. In nor¬ 
malising to the far-pair systems, positive values in Fig. 
indicate that SF is enhanced relative to the far-pair sample, 
while negative values display that it has been suppressed. We 
find that when measuring SFRs in the FIR, MIR and using 
UV-I-TIR SFRs we would see either no change or in fact 
suppression of sSFR with pair separation, while for other 
SFR indicators (FUV, NUV, O.lGyr and Ha) we would see 
net enhancement of SF. The most significant variation is 
seen in the SFRua, where sSFRs appear strongly enhanced 
by an interaction for both the intermediate- and close-pair 
samples. This echoes the findings of [Patton et al.| ( |2013| ) 
and others who find strong enhancements in emission line 
derived SFRs. 

While the results discussed above are seemingly in con¬ 
flict (or at least highlight that care must be taken when 
comparing different SFR indicators) they do not consider 
the vastly different properties of both the individual galax¬ 
ies and the merger as a whole. In the following subsections 
we shall investigate this further by splitting our samples on 
stellar mass, primary(higher mass)/secondary(lower mass) 
status, and pair mass ratio, and considering variations in SF 
as a function of pair separation directly (rather than binning 
into close, intermediate and far pairs). For completeness, in 
Appendix A we display the samples defined in [Robotham etj 
al. ( 2014[ ) but split further into major/minor mergers and 
primary/secondary status. 


4.2 Independent measures of SF as a function of 
Pair Separation, Pair Mass Ratio and Pair 
Status 

In this section, we split our pair galaxy sample into stellar 
mass bins of Alogio[M*]= 0.5 from 10®'® Mq to M©. 

We do not include galaxies outside of this mass range, where 


■ Close pairs 

■ Intermediate pairs 
□ Far pairs 


S/N 




J . i 

■D--□- 


O ■ 


SFR indicator 


Figure 5. Comparison of SFR indicators for pair galaxy samples 
discussed in [Robotham et at] | [2014[ |. Points display the median 
sSFR for each indicator, while errors show the standard error 
on the median. All points are normalised to the far-pairs value 
to highlight differences between far-, likely non-interacting, and 
close-, interacting, pairs. SFR indicators are ordered from long 
to short duration (left to right). SF appears either constant or 
increased for close pairs in comparison to far pairs depending on 
the SFR indicator used. 


sample sizes are too low to apply the analysis discussed in 
this work. We are currently investigating star-formation in 
M* < 10®'® Mq pair galaxies using a different, but comple¬ 
mentary, approach - this will be the subject a future work 
(Davies et al. in prep). 

We then bin on A rsep = 10 kpch~^ scales (i.e. ignor¬ 
ing line of sight separations which provide little information 
inside of the initial pair selections). 

As discussed previously, we calculate the excess SF for 
each indicator in comparison to a mass and redshift matched 
control sample. In order to define this control sample, for 
each pair galaxy in our sample we select three correspond¬ 
ing non-pair GAMA galaxies, which are the closest match 
in stellar mass and redshift parameter space (determined 
by the closest match in 2D space in terms of Scr clipped 
standard deviation in each parameter). We do this in each 
stellar mass, pair mass ratio, primary/secondary status and 
pair separation bin, such that each date point is scaled ap¬ 
propriately for its matched control sample. 

For each sample, mass range and radial separation bin 
we then calculate the 3cr clipped mean sSFR for both pair 
and control galaxies. We define the SFRexcesa as the ratio 
of the mean log[sSFR] in pair galaxies, divided by the mean 
log[sSFR] of non-pair galaxies, as follows: 


loglO [SFRex(M* , rsep)] 


/r(logl 0 [sSFRpairs (M* , raep)] ) 
/r(logl0 [sSFRcontrol (M* , rsep)] ) 

( 11 ) 


Note that the results discussed in this paper are not 
sensitive to choice of choice mean/median and all trends are 
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Table 2. Number of pair galaxies used in our analysis for each 
SFR indicator 


Indicator Number of pair galaxies 


FIR 

12,168 

MIR 

4,520 

UV+TIR 

10,678 

FUV 

24,285 

NUV 

24,285 

O.lGyr 

24,285 

Ha 

20,170 


observable when using median distributions over 3cr clipped 
means. We do not show the median distributions here for 
sake of clarity. Table displays the number of sources used 
in each sample. 

The left columns of Figures [6| - pd] display SFR excess for 
our full sample of pair galaxies as a function of pair sepa¬ 
ration split on stellar mass and into primary and secondary 
galaxy within the pair (green and orange lines respectively). 
The coloured polygons show the standard error on the mean 
for the pairs and non-pairs summed in quadrature at each 
radial distance. In the middle and right columns of these 
figures we split our samples further into major and minor 
mergers respectively. 

A potential pitfall of this analysis is that poor resolu¬ 
tion of continuum observations used in this work mean that 
in many cases close systems share the same resolution el¬ 
ement. While attempts have been made to deblend fluxes 
into individual sources, this may lead to systematic errors 
in our result. For instance, the secondary galaxy in a pair 
may have its emission systematically boosted/suppressed via 
deblending with the primary galaxy and vice versa. In the 
following figures we display the physical size of the PSF 
FWHM of the primary instrument used in calculating the 
star-formation rate, scaled to the median redshift of our 
pairs sample, as a dashed blue vertical line. As such, sources 
within this line may have uncertainties on their SFRs which 
arise from source confusion. However, we note that if real, 
such uncertainties in this region should be mirrored between 
primary and secondary populations, i.e. If secondary galax¬ 
ies systematically have their fluxes underestimated at close 
pair separations, then primary galaxies should systemati¬ 
cally have their fluxes enhanced, as total flux must be dis¬ 
tributed between the two sources via deblending. This effect 
is also likely to be minimal in the Magphys results which use 
the full SED to estimate SFRs. We also note that this effect 
will not be apparent in the Ha measurements, where the 
fibre aperture size is much smaller then the physical separa¬ 
tion between galaxies. In Sec. |5.2| we compare the slopes of 
these distributions. However, we note here that we consider 
slopes for both the full range of pair separations, and only 
at separations larger then the blue line (and so not affected 
by source confusion) and the results are consistent, albeit at 
lower significance. 

Below we highlight key observables from Figures [^- |11[ 
these will be discussed further in the following sections: 

• SFRfir (Fig.j^ - Firstly considering our full sample (left 
column), we find no strong enhancement or suppression of 
SFRfir as a function of pair separation - with both primary 
and secondary pair galaxies showing no real excess/deficit 


in comparison to the control sample at all pair separations. 
This is consistent with the results seen in Fig.[^ that there is 
no significant difference in FIR measure sSFR when consid¬ 
ering close-, intermediate- and far-pairs. Splitting this sam¬ 
ple into major and minor mergers, we find that the major 
merger systems (middle column) are almost identical to the 
full pairs sample and also show no excess/deficit in SF. 

However, in minor merger systems we see some dif¬ 
ferences between primary and secondary galaxies. Pri¬ 
mary galaxies show marginal enhancement of SF as a 
function of decreasing pair separation (especially in the 
10.0<logio[M*]<10.5 range - middle panel of the right col¬ 
umn) , while secondary galaxies show a deficit in SF in com¬ 
parison to the control sample. This suppression also appears 
stronger with increasing stellar mass (from top to bottom of 
the right column). We also note that, while this suppression 
is most apparent at very close pair separations, the trend of 
increasing deficit in SF as a function of pair separation be¬ 
gins outside of the PACS lOO^m PSF (blue vertical dashed 
line) and as such is not likely to be driven by deblending 
confusion. 

. SFRmir (Fig. [7]) - For all samples using our MIR SFR 
indicator we see marginal enhancement of SF as a func¬ 
tion of decreasing pair separation (all samples have a weak 
slope which increases to small pair separations). We see lit¬ 
tle difference between primary and secondary galaxies in all 
samples. 

• SFRuv+tir (Fig. - For our UV-I-TIR SFR indicator 
we find no strong enhancement/suppression of SF with de¬ 
creasing pair separation. We do see some enhancement in 
SF for major mergers in our lowest mass bin (top panel 
of the middle column), and a minor difference between 
the primary and secondary galaxies in minor mergers at 
10.0<logio[M*]<10.5 (middle panel of right column). 

• SFRnuv (Fig. 0 - For the first of our short duration 
SFR indicators (< 100 Myr) we start to see strong trends 
in our samples as a function of pair separation. For all pair 
galaxies combined (left column) we find increasing SF with 
deceasing pair separation for both primary and secondary 
galaxies. Splitting into major and minor mergers, for ma¬ 
jor mergers we find that primary and secondary galaxies 
look identical, both have increasing star-formation with de¬ 
creasing pair separation. However, for minor mergers we see 
subtle differences between primary and secondary galaxies. 
The primary galaxies show a similar relation to those in ma¬ 
jor mergers, while the secondary galaxies have no enhance¬ 
ment/suppression of SF with pair separation, except for the 
closest pairs. As with our FIR indicator, there is also a subtle 
suggestion that the normalisation of the secondary galaxies 
in minor mergers scales as a function of galaxy mass, with 
higher mass galaxies appearing more suppressed than lower 
mass galaxies (the orange line is lower in the bottom panel 
than the top panel of the right column). We do not display 
both FUV and NUV figures as both distributions are almost 
identical, but with slight normalisation scaling. 

• SFRo.iGyr (Fig. |10| ) - Our SFRo.iOyr measure displays 
similar results to the SFRnuv, but also with a slight normal- 
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Figure 6. The 3cr clipped mean excess logio[sSFRFiR] for pair galaxies as a function of pair separation binned on A rsep=10 kpch~^ 
scales. Coloured polygons show the standard error on the mean at each separation. We split out pair samples by stellar mass (rows), 
into major and minor mergers (separated at 3:1 pair mass ratio - columns) and by primary or secondary status within the pair (green 
and orange lines respectively). Blue dashed vertical line shows the FWHM of the PSF of the PACS 100/im instrument at the median 
redshift of our pairs samples. 


isation change (This is not surprising as they are all derived 
from the Magphys fits). However, for minor mergers and 
our highest mass bin (bottom right panel) we now see clear 
suppression of secondary galaxies at close pair separations. 


• SFRhq (Fig. Ill - Here we see the strongest variation 
in SFRs as a function of pair separation and the most dra¬ 
matic differences between primary and secondary galaxies 
in minor mergers. When considering all pairs (left column) 
we find an increase in SFR for both primary and secondary 
galaxies (with the enhancement much more dramatic in pri¬ 
mary galaxies). This is once again consistent with the re¬ 
sults seen in Fig. where Ha derived SFRs show a strong 
enhancement of SF when considering close-, intermediate- 


and far-pairs. However, when we consider just pair galaxies 
in minor mergers (right column), we see a dramatic change 
in the SFR with regards to primary/secondary status. Pri¬ 
mary galaxies still have their SF strongly enhanced, while 
secondary galaxies show strong suppression of SF, even at 
large pair separations (out to ~rsep >90kpc/h).This high¬ 
lights that SF is suppressed in the secondary galaxies of 
minor mergers. As SFRua measures SFRs on the shortest 
timescales, it is likely that the enhancement and suppression 
seen in primary and secondary galaxies respectively is a di¬ 
rect consequence of the galaxy interactions. We remind the 
reader that our Ha flux measurements are not subject to 
measurement errors from confusion between close pairs, and 
as such are not systematically biased by errors in flux distri- 
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Figure 7. Same as Fig. but for sSFRmir- Green dashed vertical line shows the FWHM of the PSF of W4 at the median redshift of 
our pairs samples. 


bution. Another possible explanation for the trends seen in 
Ha, are that the distribution of star-forming regions within 
a galaxy changes during the interaction. The recent simnla- 
tions of Moreno et al. (20151 show that star-formation can 
be significantly enhanced in the central regions of a galaxy 
during an interaction, but is largely suppressed in its outer 
parts. As our Ha observations only probe the central regions 
(due to the aperture-based spectroscopy) we may only be 
witnessing centrally concentrated suppression/enhancement 
of star-formation which is not a true representation of the 
global affects of the interaction. However, the trends seen in 
our Ha SFRs are consistent with those in our other short- 
duration SFR indicators, albeit at a higher significance. Fur¬ 
ther investigation into the spatial distribution of Ha SF in 
a subsample of our interacting galaxy sample is underway, 
and will be the subject of an a upcoming paper. 


• Global observation - An interesting additional global 
observation (as alluded to previously), is that pair galax¬ 
ies appear to become more suppressed in SF as a function 
of increasing stellar mass. This can be seen in all figures 
in that the lines systematically drop to lower excess sSFR 
from the top to bottom panels, and suggests that SF in 
pairs galaxies is not only affected by pair mass ratio and 
primary/secondary status, but also individual galaxy mass. 


4-2.1 Summary of Key Observables in this Section 

In summary, we find that the primary galaxies in both major 
and minor mergers show consistent trends of enhanced SF 
over our control sample as a function of decreasing pair sep¬ 
aration, and that this enhancement is more pronounced in 
the short duration SFR indicators of SFRua, SFRpuv/Nuv 
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Figure 8. Same as Fig. I^but for sSFRuv+tir- Green dashed vertical line shows the FWHM of the PSF of GALEX-NUV at the median 
redshift of our pairs samples. 


and SFRo.iGyr- In contrast, the SF in secondary galaxies 
shows different characteristics depending on pair mass ratio. 
In major mergers SF is enhanced at close pair separations, 
and follows the primary galaxies, while in minor mergers, 
we see the converse, that secondary galaxies appear to be 
suppressed at close pair separations. This effect is also more 
pronounced for short timescale SFR indictors. We find that 
at large pair separations (rgep >60 kpc/h) pair galaxies look 
similar in their SFR characteristics and as such, we are likely 
to be probing these systems prior to significant effects of 
the interaction. However, we see significant modifications to 
SFRs at close pair separations, suggesting the effect is due 
to the interaction. We also find that the normalisation of the 
distributions in Figures[^-[^marginally drops as a function 
of stellar mass, which potentially indicates that individual 
galaxy mass also has an impact on how interactions modify 
SF. 


What is clear from this analysis is that both primary 
and secondary status within the pair and pair mass ratio 
play an important role in the modification of SFRs in inter¬ 
actions and we most clearly see the effect of these parameters 
when measuring SFRs on short timescales. 


5 COMPARISONS OF SFR INDICATORS TO 
PROBE SF ON DIFFERENT TIMESCALES 

In this section we take our analysis further and combine the 
effects of both pair separation and SFR indicator timescale. 
We also use this analysis to highlight the key differences be¬ 
tween primary and secondary galaxies, and major and minor 
mergers discussed above. 
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Figure 9. Same as Fig. but for sSFRnuv- Green dashed vertical line shows the FWHM of the PSF of GALEX-NUV at the median 
redshift of our pairs samples. For our SFRpuv indicator the distributions are almost identical with a slight difference in normalisation - 
as such we do not display both figures. 


5.1 Scatter in the SFRj^a vs SFRf/h relation 

Prior to a more sophisticated analysis, and as allnded to ear¬ 
lier, Fig. 1^ highlights a possible avenne for exploring SFRs 
in close pairs using multiple SFR indicators. The SFR^a vs 
SFRf/_r relation in Fig. [^displays large scatter with signif¬ 
icant offset between Ha and FIR derived values for a large 
number of sources. As the SFRp/ij probes a much longer 
timescale than SFRpfa, could this highlight short timescale 
changes in SF for pair galaxies? 

In Fig. we show the offset between logiojSFRF/ij] 
and logio[SFR_f/a] for pair galaxies as a function of pair sepa¬ 
ration. The black lines in Fig. [^display the running median 
(solid) and upper and lower quartiles (dashed). The green 
solid line displays the Hyperfit fit to the median binned 
values, with slope, a, displayed on the plot. While only 


marginal, we do see a slight trend between logio[SFRF/i{] 
and logiojSFRrfa] as a function of pair separations. Galaxies 
at large pair separations have slightly larger logio[SFRF/i{] 
than logio [SFRpfQ,] while at small pair separations, the in¬ 
verse is true. This alludes to the fact the SF, as measured 
by Ha line emission, is very marginally enhanced over that 
measured from the FIR continuum for close pairs - poten¬ 
tially highlighting that the enhancement of SF occurs on 
timescales of < 100 Myr. However, the Spearman’s rank cor¬ 
relation coefficient for logio[SFR_Fjfl] - logic [SFRrror] against 
pair separation if just 0.0929, indicating that there is little 
statistical correlation between the two variables. We shall 
investigate this potential correlation using a more sophisti¬ 
cated method in the following section. 
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Pair Separation, kpc/h 

Figure 10. Same as Fig.j^but for sSFRo.iOyr- Green dashed vertical line shows the FWHM of the PSF of GALEX-NUV at the median 
redshift of our pairs samples - this band is the main driver of SFR calculated using Magphys. 


5.2 Direct comparison of all SFR indicators 

Taking this further, in this section we compare all of our SFR 
indicators directly and discuss the implications of variations 
between each measure. In order to make this comparison 
we use the slope of the sSFR vs pair separation relations 
in Figures[^- |ll| over the full 100kpc separation range. We 
once again use the Hyperfit package to fit the slope of 
each pair separation binned distribution and give Icr errors. 
While the distributions in Figures §-[TT] may not be well 
fit by a straight line, we stick to this simple model to avoid 
complication, and to highlight the general trend of increas¬ 
ing/decreasing excess star-formation as a function of pair 
separation. 

In Fig. we show the Hyperfit slope for all SFR in¬ 
dicators split by stellar mass, major/minor mergers and pri¬ 
mary/secondary status. The vertical blue dashed line sepa¬ 


rates long (> 100 Myr) and short (< 100 Myr) duration SFR 
indicators, as defined earlier in this work. We also highlight 
the Ha derived SFR slope in the grey shaded region, not¬ 
ing that this point is derived from nebular emission line 
measurements from aperture spectroscopy, and as such is 
subject to different biases to our continuum derived SFRs. 
The dashed horizontal line displays the dividing line between 
positive and negative slopes, this highlights either enhance¬ 
ment (slope increasing to smaller pair separations) or sup¬ 
pression (slope decreasing to smaller pair separations) of SF 
in the interaction in comparison the the most distant pairs 
(i.e. do close pairs show more or less excess SF than distant 
pairs?). As such, any significant deviation from zero slope 
displays that the interaction is affecting SF for that par¬ 
ticular sample. We note here again, that we also fit these 
distributions outside of the PSF scale for the main photo¬ 
metric band used in the analysis (blue vertical dashed lines 
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Figure 11. Same as Fig.j^but for sSFRhck- Green dashed vertical line shows the physical size of a 2 " aperture at the median redshift 
of our pairs samples. 


in Figures [^- |11| |, and find that the correlations hold true - 
but at a less significant level. 

Comparing the slopes of the distribution as a function 
of SFR indicator, we see some interesting trends. Firstly, 
considering just galaxies in major mergers, we hnd that as 
we move to SFR indicators that probe shorter timescales, 
star-formation is enhanced more strongly (for long duration 
measures there is no star-formation enhancement, while for 
short duration measures we see consistent enhancement of 
star-formation). This is true for both primary and secondary 
galaxies. It is also interesting to note that once again there 
appears to be a weak trend with stellar mass, however, this is 
only apparent in secondary galaxies, with secondary galax¬ 
ies being less enhanced at increasing stellar mass (the or¬ 
ange points drop from the top to bottom panels). In con¬ 
trast for galaxies in minor mergers, we find that primary 
and secondary galaxies show different characteristics. Pri¬ 


mary galaxies display similar trends to the primary galaxies 
in major mergers, showing excess star-formation for short 
duration measures. Secondary galaxies largely show suppres¬ 
sion of star-formation, in all but the smallest mass bin. What 
is clear from this hgure is that for short duration SFR indi¬ 
cators, secondary galaxies in minor mergers show different 
characteristics to their primary counterparts. 

These results are indicative of SF being enhanced in pri¬ 
mary galaxies of all mergers and secondary galaxies being 
enhanced in major mergers and suppressed in minor merg¬ 
ers, but only when measured on short timescales. 

To highlight the dependance on primary/secondary sta¬ 
tus and pair mass ratio of the modification of SFRs in inter¬ 
actions further, we compare differences between the slopes 
for major and minor mergers, at fixed stellar mass binning 
and pair status (Fig. |14[ | and differences between primary 
and secondary galaxies at fixed pair mass ratio (Fig. |15[). 
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galaxies appear suppressed relative to primaries for all star- 
formation rate indicators and in the majority of stellar mass 
bins (the majority of points in the right panel of Fig. 15 sit 
above the central line). This effect also appears more pro¬ 
nounced for short duration indicators and at intermediate 
stellar masses. 

In combination, these figures highlight that the sec¬ 
ondary galaxies in minor mergers show distinctly different 
SF characteristics to all other pair galaxies, but this ef¬ 
fect is only strong when SFRs are measured on < 100 Myr 
timescales. 


6 SUMMARY OF THE EFFECT OF CLOSE 
INTERACTIONS ON GALAXY SFRS AND 
DISCUSSION 


Pair Separation, Kpc 


Figure 12. logiojSFR^/jj] - logio[SFR//(j] for our pairs sample 
as a function of pair separation. Black lines show the running me¬ 
dian (solid) with upper and lower quartiles (dashed). The green 
solid line displays the Hyperfit fit to the median binned distri¬ 
bution (with slope, a). While marginal we see a slight trend in 
logic[SFRi^’/ij] - logic [SFRuor] suggesting when considering all 
pair galaxies together, SF may be enhanced on short timescales 
in closely interacting pairs. 


Error bars in these figures are derived from the sum of the 
squares of the Hyperfit errors given in Fig. In Fig. [Td} 
offsets from the central line display a significant difference in 
the effect of interactions between major and minor mergers, 
where positive values show that SF is more strongly sup¬ 
pressed in a minor merger. Note that this figure does not 
show absolute values and as such zero does not highlight 
that there is no change to SFRs in the interaction, only that 
there is no difference between major/minor mergers galax¬ 
ies. We find that primary galaxies are largely consistent with 
their being no difference between major and minor mergers 
for all SFR indicators (the green points are largely consis¬ 
tent with the central line given the errors), while the sec¬ 
ondary galaxies stronger suppression of SF in minor mergers 
than major mergers (the orange points sit above the central 
line). This effect appears largest in short duration SFR in¬ 
dicators (i.e. the differences between secondary galaxies in 
major and minor mergers is most apparent when probed on 
short timescales). 

In Fig. [^offsets from the central line display a signifi¬ 
cant difference in the effect of interactions between primary 
and secondary galaxies, where positive values show that SF 
is more strongly suppressed in the secondary galaxy. We 
Hnd that for major mergers there is little difference between 
primary and secondary galaxies specifically for long dura¬ 
tion indicators (as discussed previously). We do see that 
for short duration timescale indicators, secondary galaxies 
are marginally suppressed in comparison to primary galax¬ 
ies when considering the higher stellar mass bins only (the 
gold and purple triangles sit above the central line), but this 
effect is small. For minor mergers we see a larger discrep¬ 
ancy between primary and secondary galaxies. Secondary 


Now we outline the key results outlined in the previous two 
sections and discuss a possible toy model for SF in galaxy 
mergers which may explain these results: 


• Considering pair populations as a whole, SF can appear 
either unaffected or enhanced by an interaction depending 
on SFR indicator used (Fig. [^. 

• At large pair separations (>50kpc) we see little differ¬ 
ence between different galaxy populations, suggesting these 
systems are not being strongly effected by the interaction 
(Figures 

• At close pair separations (<30kpc) we see both enhance¬ 
ment and suppression of galaxy SFRs depending on the sub¬ 
sample of galaxies we select, and this modification to SFRs 
appears more strongly for short duration SFR indicators 
(Figures 


• For long timescale SFR indicators, we see little enhance¬ 
ment or suppression of SF in interacting galaxies. This holds 
true for both major and minor mergers, and for both pri¬ 
mary and secondary galaxies - suggesting that we are prob¬ 
ing the SF in the galaxy prior to its SFR being modihed by 


the interaction (Fig. 131. 


• For short timescale SFR indicators (especially Ha), pri¬ 
mary galaxies show enhancement in SF at close pair sepa¬ 
rations, which is consistent across major and minor mergers 

(Fig. ft. 


• For short timescale SFR indicators, secondary galaxies 
show distinctly different characteristics between major and 
minor mergers. In major mergers SF is enhanced with de¬ 
creasing pair separation, while in minor mergers SF is sup¬ 


pressed with decreasing pair separation (Figures 13 14 & 


• Primary galaxies appear largely agnostic to pair mass ratio 
while secondary galaxies show suppression in minor mergers 


only (Figures 13 14 & 151. 


To relate these observables to the true effects of interac¬ 
tions on SFRs, we use both the pair separation and variation 
of SFR indicators on different timescales as a proxy for the 
stage of the galaxy interaction. We relate both our long du¬ 
ration SFR indicators at all separations and our short dura¬ 
tion measures at large pair separations to be representative 
of the pre-, or early-interaction stage. Conversely, our short 
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Figure 13. Comparison of the slopes of logio[sSFR] vs pair separation for all SFR indicators shown in Figure [^- |11| Samples are split 
on stellar mass (rows) and pair mass ratio (columns), as well as primary (green) and secondary (orange) status within the pair. We 
divide of SFR indicators into long and short duration with the vertical blue dashed line, and highlight the Hq emission line derived SFR 
with the grey shaded region. Significant deviations from the central horizontal line display modifications to SF induced by the a merger 
- positive values display enhancement of SF and negative values suppression of SF. Error bars display the la error on the hyper.fit fit to 
the slope. 


duration measures at close pair separations is representative 
of the late-interaction stage. This simplistic model assumes 
that galaxies are all on their first approach and that they 
have not had previous passes, which would effect their star- 
formation. Such systems could reach relatively large pair 
separations, but would have already had their star-formation 
enhanced/suppressed by the interaction. In our analysis we 
have no method of determining if galaxies have already had 
one or more passes prior to the observation epoch, as such 
we can not account for this affect. However, we note that 


Robotham et al. (20141 find that, for the pair sample used 


in this work, galaxies at r^ep < 20kpc show much greater 
signs of visual disturbance than those at 20 < r^ep < 50 kpc 
and both these populations show higher visual disturbance 
than pairs at rsep > 50 kpc. This indicates that galaxies 
are more strongly affected by interactions at closer separa¬ 
tions, and while a fraction of galaxies at large separations 


my have already had a number of passes, the major trend is 
that closer separations mean stronger pair interactions and 
a later stage in the merger process. Therefore, we deem it is 
reasonable to assume that to first order pair separation can 
be directly related to merger timescale. 

Using these assumptions we can piece together a hy¬ 
pothetical toy model for SF in the primary and secondary 
galaxies in both major and minor merger systems. Fig. 
displays this hypothetical model with interaction time/pair 
separation plotted against current SFR. We also display the 
timescales over which both long and short duration SFRs are 
measured, for two different scenarios: tl, when the galaxies 
are at large pair separations (an early stage in their interac¬ 
tion) and t2, when the galaxies are at small pair separations 
(a late stage in their interaction). 

In our major merger model SF remains relatively con¬ 
stant throughout the merger in both the primary galaxy and 
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Figure 14. Comparison of the difference of the logio[sSFR] vs pair separation slope between major and minor mergers. Values given are 
the offset between corresponding points in the left and right column in Fig. |13| Significant offsets from the central line show differences 
between major and minor mergers for each population and SFR indicator. Once again, secondary galaxies in show more significant 
differences between major and minor mergers. 
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Figure 15. Comparison of the difference of the logic [sSFR] vs pair separation slope between primary and secondary galaxies. Values 
given are the offset between primary and secondary galaxies in each row in Fig. |13| Significant offsets from the central line show differences 
between primary and secondary galaxies at each stellar mass and SFR indicator. Once again, secondary galaxies in minor mergers in 
show more significant differences between primary and secondary galaxies, specifically for short duration SFR indicators. 
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secondary galaxy, but is enhanced at the very late stages (we 
see the strongest enhancement in the shortest timescale in¬ 
dicator - Ha). These enhancements in SF are likely to be be 
due to gas turbulence, tidal torques or shocks cause by the 
gravitational attraction of the other system. In this scenario, 
the secondary galaxy has sufficient mass (in comparison to 
the primary) to retain its gas and continue star-forming at 
a relatively constant rate i.e. it is not tidally stripped. This 
is consistent with our observations summarised above. The 
primary and secondary galaxies show some enhancement of 
SF with pair separation, and this increases for the short¬ 
est timescale indicator. If SFRs were measured at tl and t2 
we would observe a marginally increasing measure of both 
long and short duration SFRs (the primary and secondary 
galaxy show similar SF properties which increase in SF for 
short duration measures). The largest tidal interaction will 
occur at the smallest pair separation and be most evident 
in the shortest duration SFR indicator, as such we see the 
largest enhancement of SF in Ha derived SFRs. 


For our minor merger model, the primary galaxy fol¬ 
lows the same SF evolution as both systems in the major 
merger. It is sufficiently high mass in comparison to the sec¬ 
ondary to retain its gas and continue star-forming, and has 
its SFR enhanced over time as interactions with the sec¬ 
ondary systems cause gas turbulence and shocks. However, 
the secondary system follows a different SF evolution path 
during the interaction. Initially its SF is not significantly 
altered and hence we do not see distinct differences in long 
duration SFR measures or a large pair separations, which 
probe this epoch. However, at some stage during the interac¬ 
tion SF begins to be suppressed. This is likely through tidal 
stripping or gas heating. This is evident as the suppressed 
SFRs at close pair separations and in short duration SFR 
indicators. If we were to measure the SFR of the secondary 
system at tl and t2, we would find that at tl (large sepa¬ 
ration), the galaxy would look somewhat like the primary 
galaxy, and we would see no suppression of SF. However, at 
t2 we would find short duration SFRs suppressed relative 
to tl, this effect would not yet be observable in the long 
duration SFR measures. At this point the secondary galaxy 
SFRs would be distinctly different from the primary galaxy 
when observed using short duration indicators - this model 
is also consistent with the observations summarised above. 


7 CONCLUSIONS 

We have investigated star-formation activity as a function 
of pair separation for interacting galaxies for different 
merger scenarios, using SFRs derived from multiple SFR 
indicators covering a range of timescales. The key results 
from our analysis are as follows: 

• SF in interactions proceeds differently for galaxies in 
minor mergers than in major mergers. 

• SF in minor mergers proceeds differently in galaxies 
depending on whether they are the primary (high mass) or 
secondary (low mass) galaxy. 

• Primary galaxies show some enhancement of SF in 
all mergers which is only evident the shortest timescales - 
and as such the later stages of the interaction. 

• Secondary galaxies show little change in SF in ma¬ 
jor mergers, except for some enhancement on the shortest 
timescales, but are suppressed in minor mergers. 

• SF suppression in secondary galaxies occurs on timescales 
of < 100 Myr as it is not evident in long duration SFR 
measures - as such, it also occurs at a late stage in the 
interaction. 


Using these results we propose a scenario for the evolu¬ 
tion of SF in major and minor merger systems. In the major 
mergers both primary and secondary galaxies have their SF 
marginally enhanced through tidal turbulence and shocks 
for the early stages. At the later stages of the merger SF is 
likely to be enhanced. SF is not suppressed as both galaxies 
have sufficient mass (in relation to the other) to retain their 
gas. In the minor mergers the primary galaxy follows the 
same SF enhancement as the major merger systems, while 
the secondary galaxy has its SF suppressed as its gas is ei¬ 
ther tidally stripped, heated or simply stretched to a lower 
mean density which is not sufficient for SF. This toy model 
is consistent with the recent findings of [Robotham et ah] 
(20131 and De Propris et al. (20141, as well as the galaxy 
merger simulation of |Di Matteo et al.| ( |2007[ ) - predicting 
that the effects of galaxy interactions on star-formation are 
dependant on both pair mass ratio and primary/secondary 
status within the pair. 


If representative of the true galaxy-galaxy interactions, 
these toy models could potentially explain the contradic¬ 
tory results from previous studies which find both SF en¬ 
hancement and suppression in close pair systems. Our mod¬ 
els predict that SF is enhanced in the majority of cases - 
where mass ratios are close to 1:1, as well as in the primary 
galaxy in a minor merger. However, SF is suppressed in the 
secondary galaxy in minor mergers, but only at close pair 
separations. As such, whether SF is enhanced or suppressed 
in galaxy-galaxy interactions is largely dependent on the 
pair mass ratio and pair separation. Clearly by observing 
galaxies at different interaction stages and at different pair 
mass ratios, we will will obtain conflicting results as to the 
suppression/enhancement of SF in interacting systems. 
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APPENDIX A: ROBOTHAM ET AL SAMPLE 
SPLIT ON MAJOR/MINOR MERGER AND 
PRIMARY/SECONDARY STATUS 

In Fig. |A1| we display the pair sample of [Robotham et al.| 
| |2014| ) discussed in Sec.|^but split into major/minor merg¬ 
ers and primary/secondary status. These distributions echo 
the main results derived in this paper. Gonsidering the pop¬ 
ulations as a whole, we see enhancement of star-format ion in 
close pairs, but this enhancement is only observable in SFR 
indicators which probe short timescales. Splitting the sam¬ 
ples into primary and secondary galaxies, we find that the 
increase in SF for close pairs is completely driven by the pri¬ 
mary galaxies, where strong excess is seen for short duration 
indicators (specifically Ha). Splitting further on both pair 
mass ration and galaxy status, we find that it is the primary 
galaxies in both major and minor mergers that show strong 
enhancement in SF. For secondary galaxies we see little dif¬ 
ference between far-, intermediate- and close-pairs for all 
SFR indicators. This is consistent with the main results in 
this paper, when derived in a separate, but complementary, 
manner. 
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Figure Al. Pair galaxies selected from [Robotham et al.| ( |2014| | split by major/minor mergers and primary/secondary status. All pair 
samples are normalised to the far-pairs’ value to highlight the effects of interactions on pair galaxies in each SFR indicator. 













































